A Dynamic Model of Streamer Coupling for High Pressure Discharges 
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A streamer coupling theory is developed to describe the formation of homogenous emission, and 
the high moving speed of emission patterns in high pressure discharges. By considering the effects of 
both electron diffusion and electronic drift in the streamer head, the minimum required preionization 
level n m in for the formation of streamer coupling is found to depend on electric field strength, gas 
pressure and electron temperature. The homogeneity and moving speed of the emission pattern in 
streamer coupling head increase with preionization level no, when no > n m i n . The predicted results 



for atmospheric helium plasma indicate n n 
agreement with experiments. 



10 cm and moving speed of 10 — 10 m/s, in 



PACS numbers: 51.50.+V, 52.80.Dy, 52.80.Hc, 52.80.Mg 



Homogeneous or glow-like emission in high pressure 
discharges has been studied for decades, and most of the 
works focused on experimental conditions for their for- 
mation and methods of improving discharge homogene- 
ity [lfl ■ Although high preionization level is verified to be 
necessary for homogeneous discharges [lJ-[3] , its forma- 
tion mechanism, different from those of glow discharges 
in low pressure and streamer discharges in high pressure, 
has not been fully clarified [3|, |j|. In a recent experi- 
ment, atmospheric pressure plasma jet (APPJ), originat- 
ing from dielectric barrier discharge (DBD) and spray- 
ing into ambient air, is found to be a traveling bullctlikc 
plasma volume with a high moving speed of the order of 
10 4 ~ 10 6 m/s and a ring-shaped cross-sectional emission 
pattern, named as "plasma bullet" [^-B- A model of self- 
sustained photo-ionization streamer [fj|, first developed 
by Dawson and Winn [l(| , was invoked to explain the na- 
ture of plasma bullet as one single streamer [lfj, L!J| • Al- 
though the model sufficed to account for the scale of mov- 
ing speed, it was inadequate to explain the ring-shaped 
emission pattern of axial-symmetrical homogeneity and 
the change of moving speed with external electric field 
0-Q. Thus, we suppose that APPJ is neither a sin- 
gle streamer discharge, nor traditional Townsend or glow 
discharges in low pressure, which are homogeneous in ra- 
dial direction and brightest in the vicinity of the anode 
or cathode [ll-y|.The discharge mechanism should be a 
new one, which we here call it streamer coupling. 

The thoughts of streamer coupling model was primar- 
ily advocated by Palmer to predict a volume-stabilized 
glow-like discharge in atmospheric pressure helium dis- 
charge 12]. In Palmer's theory, the interaction of simul- 
taneous developing streamers leads to the formation of 
one large discharge canal, and the dominating force in 
each streamer head is electron diffusion yjj [l4(. How- 
ever, experimental results show that, in high pressure 
experiments such as atmospheric pressure discharge, the 
dominating force responsible for the electron cloud ex- 
pansion in a streamer head is the electrostatic repulsion 
of high-density charged particles, instead of the diffu- 
sion caused by electron density gradient (l5|. Therefore 



in Palmer's model the predicted minimum preionization 
level of 10 4 cm" 3 for homogeneous discharge [12[ is not 
consistent with the experimental value of 10 5 cm -3 [14j . 
Another defect in Palmer's model is that there is no ex- 
plicit relationship between important physical properties 
of the discharge gas with the formation of homogeneous 
discharge, such as electric field and gas pressure. In this 
work, by considering the electron diffusion, electrostatic 
repulsion in streamer head, and electron drift under the 
electric field, an improved streamer coupling model is 
proposed to describe the dynamics of high pressure dis- 
charge patterns. 

We consider the "fluid approximation" for each 
streamer. For the sake of simplicity, we here investi- 
gate only the primary anode-directed streamer and as- 
sume it propagates in a uniform background electric field 
Eq . The continuity equation for the basic dynamics of a 
streamer formation and propagation is, 



-^r = ~ V ' ( n eVe) + S e , 
at 



(1) 



where n e is the electron density, S e is the electron source 
term, and v e is the electron velocity, determined by 
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with fi e and D e are the electron mobility and the diffusion 
coefficient, respectively. A detailed analysis responsible 
for these equations is given as follows. 

(1) Resulting from drift and diffusion of charged par- 
ticles in the local electric field E, streamer propagation 
is mainly determined by the motion of electrons. The 
ions can be treated as immovable particles, since their 
mobility \Xi and diffusion coefficient Di can actually be 
neglected, comparing with those of electrons jla - tl7l |. 

(2) The source term S e can be treated in two processes. 
In the process of primary avalanche, the source term S e is 
directly proportional to the exponent of first Townsend 
ionization rate a , as S e oc e ax , where x is the length 
of the avalanche [15|. In the other process of streamer 




FIG. 1: Simultaneous avalanche-streamer transition. (a) 
Avalanches start to develop when electric field arises to a 
certain value Eq. (b) Primary avalanches turn into streamers 
when they satisfy the criterion of streamer formation at the 
length of Xq. (c) Adjacent streamers just overlap with each 
other to form streamer coupling. 



propagation, the source term S e could be treated as S e oc 
Ee~\ E °' E ' , where Eq and E are external and total electric 
fields, respectively [ll|. In the present work, we focus on 
the former case. 

(3) We concentrate on the streamer dynamics under 
the strong external electric field Eq, as in high pressure 
discharges. The criterion of streamer formation says, a 
streamer is born of an avalanche if the electric field E' in- 
duced from the space charge in the streamer head reaches 
the order of external field Eq [15| . The correspondingly 
approximate equality is, 
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Eq, 



(3) 



where Rq is the characteristic radius of space charge 
in the streamer head at the transformation point. The 
streamer head region of intensive ionization, moving to- 
gether with a strong field E = Eq + E', transforms the 
gas to plasma. A plasma channel is left due to the pro- 
duction of new plasma region. 

For the case of high pressure discharges in capacitively 
coupled plasmas, simultaneous electrons leave the dielec- 
tric surface coated on the instantaneous cathode towards 
the anode 2, 3] , when the electrode polarity connected to 
the squared power source turns from positive to negative 
half-cycle at time t = —At. At is assumed to be the aris- 
ing time of external electric field from to the value at 
which the primary avalanche starts at the leaving place of 
x = Ax, the distance to dielectric surface. After the time 
t = 0, the supplying pulsed power source is assumed to 
be sustained at constant value Eq for the streamer prop- 
agation [see Fig. Ufa)]. At a certain place x = xq with 
t = to, these simultaneous primary avalanches transform 
into simultaneous streamers. The streamer head radii in 



propagating and transverse directions are assumed to be 
R p and R t , respectively [see Figs. QJb) and (c)]. These 
simultaneous streamers will overlap if the transverse ra- 
dius Rt is larger than the half distance between adjacent 



streamer head centers (47rno/3)' 
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i.e., 
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(4) 



For better discussion bellow, the surface density of simul- 
taneous primary electrons is replaced by volume density 
no, which is known as preionization level. The equal sign 
case of Eq. ((4]) is shown in Fig. [He). The overlapping 
streamers hereafter is called streamer coupling, whose 
propagation likes one single streamer, except that it has 
a much larger space charge head. 

For the avalanche-streamer transition, if the expansive 
force at streamer head edge is dominated by electrostatic 
repulsion or electron diffusion, according to Eq. ([2]), the 
speed of electron can be expressed as: 
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In experimental breakdown condition of high pressure 
discharges, the expansion of an avalanche head is mainly 
due to repulsive force rather than diffusion one. The 
difference of the two forces in magnitude can be one or 
two orders in many cases, such as in atmospheric pressure 
air discharge fl5j . 

Due to the cancelation of induced repulsive field be- 
tween the simultaneously developed adjacent streamer 
heads in transverse direction, the dominator of Eq. ([5]) 
for the directions of propagation and transverse are re- 
pulsion and diffusion, respectively. Using Eq. ([3]) to the 
propagating direction along the avalanche development, 
we obtain: 
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(6) 



For the sake of simplicity, the assumption a — a(E ) is 
used in Eq. ((6]) during the development of an avalanche 
when the external field is only sightly distorted [15]. Also 
the corresponding transverse radius Rt increased by dif- 
fusion is: 
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where E w Eq is used from x = to x = xq in Eqs. ([6]) 
and ([7]) in the avalanche development. A more strict cal- 
culation should consider the integral of E. However, con- 
sidering the uncertainty due to no clear plasma edge like 
solid, the above approximation is enough for our estima- 
tion. Using a dimensionless streamer density £t = SUf tiq 
to the transverse direction, we have: 
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FIG. 2: Distribution of minimum required seed electron den- 
sity n m in for atmospheric pressure helium plasma with elec- 
tron temperature T e and electric field Eq. 



The criterion for streamer coupling formation can be re- 
expressed as: 



6>1, 



(9) 



and thus we obtain the expression for primary electron 

density, 
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(10) 



where Einstein relation of D e / u e = ksT e /e is used, ks is 
the Boltzmann constant, and T e is the electron temper- 
ature. n m ; n is defined as the minimum required preion- 
ization level. Since a — a(£'o,P), "-mm is a function of 
the external electric field Eq, the gas pressure P and the 
electron temperature T e . 

For the propagation of streamer coupling, the electric 
field strength at the front of streamer head Em is esti- 
mated as: 



E M 



l + ^n Q jEo if£ P >l, 
1.5Eo ' if i P < 1, 



(11) 



where £ p = ^-R^uq. The strength of electric field Em 
results from the total effect of streamer heads when £ p > 
1, and this effect disappears when £ p < 1. 

Applying the above theoretical results to the atmo- 
spheric pressure helium plasma, we find that the mini- 
mum required seed electron density is relatively low for 
the discharge of a high electron temperature and a low 
electric field (see Fig. [2]). Choosing data a w 5.3 x 
10 3 m _1 from the estimation of experimental value (18J . 
and a typical experimental condition of T e = 2 eV, Eq — 
4 kV/cm, we can obtain: 

n m f n - 1.1 x 10 5 cm" 3 

The above result is about one order of magnitude hi ghe r 
than the predicted result from Palmer's model |12l |. 
While the experimental minimum required density for 
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FIG. 3: Distribution of total electric field E projected in 
streamer direction at the front of streamers for different pri- 
mary seed electron density. We choose electron temperature 
of 2 eV and helium plasma of atmospheric pressure. 



homogeneous discharge is the order of 10 5 cm 3 |14l. |20|. 
which is in favor our predicted result. 

Based on the above calculation, the predicted result for 
streamer head radius in propagating direction is R^ c ~ 
0.03 cm, and in transverse direction is i?^° ~ 0.01 cm 
which is the radius of a single streamer. The two values 
suggest that the distribution of space charge in streamer 
head likes a " goose egg" , as shown in Fig. [TJc) . 

The transverse distribution of the total field strength 
E in the propagating direction at the front of streamer 
heads is shown in Fig. [3] Streamers do not overlap with 
each other when the primary seed electron density Uq 
is lower than the minimum required preionization level 
"-min, such as n = 0.08, 0.5 and 0.8 x 10 5 cm~ 3 . The 
electric field is continuous for the streamer coupling when 
no > n m i n , and its relative smoothness increases with no, 
such as no = n m j n , 5n m j n and 20n m j n . This indicates 
that the distribution of ionization and radiative processes 
are almost homogeneous for the streamer coupling, and 
the emission homogeneity is improved by increasing the 
preionization level no- This prediction is qualitatively 
consistent with the experimental results [3j, Il4j . which 
suggest that the homogeneous discharge can be only ob- 
tained with high preionization level. 

The development of the streamer coupling is led by the 
drift of electrons at the front of the streamer coupling 
head, since streamers propagate along the direction of 
the strongest electric field [15j . Therefore we can obtain 
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FIG. 4: Distribution of the moving speed of emission pat- 
tern with primary electron density at T e — 2 eV and Eq — 
4 kV/cm for atmospheric pressure plasma. The left and right 
regions separated by vertical dashed lines denote discharges 
of single streamers and streamer coupling, respectively. 



the moving speed of the streamer coupling head: 



f fi e E (l + ^n )E if&>l, 
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(12) 



Since the streamer coupling head is the most intensive 
ionization region, the propagation of the streamer cou- 
pling head represents the moving of discharge pattern in 
actual experiment. Using the data of electron mobility 
He = 1.1 x 10 3 cin 2 V _1 s _1 in the atmospheric pressure 
helium plasma [4j, we obtain the moving speed of the 
discharge pattern which is shown in Fig. 2] Different 
discharge regions are separated by vertical dashed lines 
according to preionization level. The left and right re- 
gions are single streamers, where streamers are almost 
independent, and the streamer coupling, where stream- 
ers overlap with each other, respectively. In the middle 
region, although the discharge is also separated stream- 
ers, the effect of other ones in transverse direction can 
not be ignored. It shows that the moving speed of the 
discharge pattern increases linearly with no and Eq in the 



streamer coupling region. The scale of the moving speed 
is consistent with the experimental results of "plasma 
bullet" 0-11]. Furthermore, according to Eq. (fT2"]l. the 
ionization rate a and primary seed electron density no 
is axisymmetrically distributed due to the axisymmet- 
rically distributed of gas mixing and dielectric surface. 
The ionization rate a reaches its maximum value when 
the content ratio of nitrogen in the helium plasma is at 
the level of 10 -3 0, 0j|. Therefore, the ring-shaped pat- 
tern of plasma bullet can also be explained by Figs. [3] 
and[U 

Although the predicted results are consistent with done 
experiments, an identifying experiment needs to be done 
to check the predictions quantitatively. For the atmo- 
spheric pressure helium discharge, by using the trans- 
versely excited atmospheric pressure CO2 laser system 
|2Cf , we can control the preionization level no in the setup 
of capacitively coupled plasma to identify Figs. |3Jand|l] 
Fig. [3] suggests that the contrast of streamer emissions 
and background increases with no when no < n m ; n , and 
the homogeneity and intensity of emission pattern in- 
creases with no when no > n m - m . Fig. 2] suggests the 
moving speed of discharge pattern under certain electric 
field and preionition level, and that it increases with no- 

To conclude, we have derived analytically the mov- 
ing speed of emission pattern and the minimum required 
preionization level for anode-directed streamer coupling, 
and supposed that the streamer coupling is required to 
generate a homogeneous emission pattern in high pres- 
sure discharges. Both values depend on the electric field 
strength and the gas pressure and the electron tempera- 
ture. Based on these predictions, we investigate the emis- 
sion homogeneity and its moving speed in atmospheric 
pressure helium plasma. Our predictions are consistent 
with experimental results. The model of the streamer 
coupling is very useful for understanding the dynamical 
process of high pressure plasma. 
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